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' Abstract 

/\ I From an analysis of the ionisation energy loss of charged particles selected from 

c3 ■ 110326 e~^e~—^T~^T~ candidates recorded by the OPAL detector at e"'"e~ centre-of- 

mass energies near the Z° resonance, we determine the one-prong tau decay branching 
ratios: 

Br(r-^z/^K- > Oh°) = 1.528 ± 0.039 ± 0.040% 

Br(r-^z/^K-) = 0.658 ± 0.024 ± 0.029% 

where the h° notation refers to a 7r°, an t], a Kg, or a Kl, and where the first uncertainty 
is statistical and the second is systematic. 
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1 Introduction 



Precise studies of the decays of the tau lepton have been made possible by the 
availability of large, low background tau-pair samples, such as those produced in e"'"e~ 
collisions at CESR and at LEP I. This paper reports on an analysis of one-prong tau 
decay modes containing a charged kaon, using the complete set of data collected by the 
OPAL experiment between 1990 and 1995 at e"'"e~ centre-of-mass energies near the Z° 
resonance. The excellent charged particle identification capability of the OPAL detector 



is exploited to obtain precise measurements of the r~— ;>//,- K~ and t^^z/t-K" > Oh° 
branching ratios, where the charge conjugate decays are imphed in these interactions 
and throughout this paper. The h° notation refers to a 7r°, an rj, or a K°, where the 
K° can be either a Kg or a Kl. 

Measurements of the r^^ z/^ K~ branching fraction may be used to determine the 
charged kaon decay constant, /k, and the CKM matrix element, {Vusl- In addition, 
precise measurements of this branching ratio may also be used to set limits on physics 
beyond the Standard Model, such as charged Higgs effects in tau decays M. 



2 The OPAL Detector 



The OPAL detector is described in detail in reference 0. The component of the 
detector most important to this analysis is a large volume jet chamber, which measures 
the momentum and specific ionisation energy loss, dE/dx, of charged particles. The jet 
chamber provides the only means of distinguishing between charged pions and kaons 
in the momentum range of interest to this study. 

The jet chamber is a cylinder 4 m long and 3.7 m in diameter, and is divided by 
cathode wire planes into 24 azimuthal sectors. Each sector contains one radial plane 
of anode wires, which are staggered to resolve left-right ambiguities. The chamber is 
contained in a solenoidal magnetic field of 0.435 T, and is filled with an argon-methane- 
isobutane gas mixture at a pressure of 4 atmospheres. This arrangement provides an 
intrinsic transverse spatial hit resolution of axy = 130 /xm, and a two-hit resolution of 
2.5 mm.Q In the barrel region of the jet chamber (| cos^| < 0.68), the ionisation energy 
loss of a charged particle is sampled up to 159 times. A truncated mean is formed by 
discarding the highest 30% of the measurements, resulting in a dE/dx resolution of 
3.1% for isolated tracks in the chamber [Q]. 

Other tracking detectors include a high precision silicon microvertex detector sur- 
rounding the beam pipe, a precision gas vertex detector, and a layer of wire chambers 
with drift direction parallel to the z-axis lying immediately outside the jet chamber. 
These ^-chambers extend to | cos^| = 0.80, which encompasses the acceptance of this 
analysis, and accurately determine the polar angle of charged particles traversing the 
central detector. A lead-glass electromagnetic calorimeter and presampler chambers 
are located outside the magnetic coil and jet chamber pressure vessel. The electromag- 
netic calorimeter is primarily used in this analysis to identify tau decays that include 
a 7r° or a K° in the final state. The return yoke of the OPAL magnet is instrumented 



^ A spherical coordinate system is used, with the +z-axis in the direction of the circulating electron 
beam. The angle 9 is defined as the polar angle with respect to the +z-axis, and <j) is defined as the 
azimuthal angle measured from the +x-axis, which points towards the centre of the LEP ring. 



for hadron calorimetery, and is surrounded by external muon chambers. The hadron 
calorimeter and muon chambers are mainly used in this analysis to select control sam- 
ples containing muons and to veto non-tau backgrounds in the tau-pair candidate 
sample. 



2.1 Parameterisation of lonisation Energy Loss 

The mean ionisation energy lost per unit distance by a particle of charge q and speed 
(3 traversing a gas is approximately described by the Bethe-Bloch equation [Q ^ : 

I^pred = -^ [A2 + In (7'/?') - /?' - 5(/?)/2] , (1) 

where Ai and A2 are constants which depend on the gas composition, and 5{13) is a 



function which describes polarisation effects in the gas. 

The ionisation energy loss parameterisation used by OPAL, -Dpred, is a variation of 
the Sternheimer-Peierls modification to the Bethe-Bloch formula p, ^, and uses six 
input parameters. Figure |I|(a) shows the dependence of dE/dx on the momentum of 
various particles in the OPAL jet chamber, as predicted by -Dpred- Figure |l|(b) shows 
the particle resolving power, TZij, versus momentum for various pairs of particle species 
i and j, where: 

/v-ij = l-^pred ~ -^prcdl / ^J 

and where a is the resolution of the measured dE'/dx. Note that the OPAL jet chamber 
yields a pion/kaon separation of at least 2a for particles between about 2 to 30 GeV/c. 

The parameters used in the calculation of -Dpred and a are tuned to the measured 
energy loss of charged particles in hadronic 7P — »• qq decays, yielding a dE/dx param- 
eterisation that is accurate enough for nearly all analyses of OPAL data. 

However, hadronic Z° -^ qq decays tend to have a much higher track multiplicity 
than one-prong tau decays, leading to small systematic differences between the mea- 
sured dE/dx in the two environments. Although these differences are less than one a 
in magnitude, they must be corrected for accurate particle identification in one-prong 
tau decays. Based on a study of the measured energy loss of muons in e+e~— >;U+/i~, 
T~ -^ fi~T7 ^ i>T- , and e^e~— >e^e~/i^/i~ control samples, the following corrections are ap- 
plied to the the OPAL dE/dx parameterisation, Dpred, and to the measured energy 
loss and its uncertainty, -Drncas and a: 

Multiplicative correction to -Dpred '■ To obtain a parameterisation which correctly pre- 
dicts the measured energy loss of particles from tau decays, -Dpred must be cor- 
rected with a multiplicative factor, s{f3). Parameterising s{(3) as a second order 



polynomial in — log(l — (3"^) yields satisfactory results. The magnitude of the 
correction is of order 1% for tracks in one-prong tau decays. 

Additive correction to -Dmeas : The measured energy loss of charged particles in low- 
multiphcity events is found to depend strongly on the azimuthal separation, 0, 
between the track and anode plane in the jet chamber cell. To improve the 
dE/dx resolution, this behaviour is corrected with a function, /(0), with seven 
parameters tuned to the measured energy loss of e+e~— »• /i^/i~ candidates. The 
magnitude of this correction is also of order 1%. 

Multiplicative correction to a: To correctly predict the spread in the dE/dx of charged 
particles from one-prong tau decays, the OPAL parameterisation of a must be 
corrected with a multiplicative factor, Sj-es ~ 0.9. 



3 Monte Carlo Generated Event Samples 



For this analysis, tau lepton decays are simulated with the KORALZ 4.0 Monte 
Carlo generator and the TAUOLA 2.4 decay package [0, p. A Monte Carlo sample of 
600000 tau-pair events is generated with input branching ratios determined from world 
averages or theoretical expectations [^. Background contributions from non-r sources 
were evaluated using Monte Carlo samples obtained from the following generators: 
hadronic Z° — > qq events were simulated using JETSET 7.4 [|I^, e"'"e~— i>/i"'"/i~ events 



using KORALZ, Bhabha events using BABAMC |11| and BHWIDE [12 , and four- 



fermion events using VERMASEREN 1.01 and FERMISV |4|. 



All Monte Carlo samples are passed through a detailed simulation of the OPAL 



detector ||15[ and are subjected to the same analysis chain as the data. 



4 Data Selection 



This analysis uses the full data set collected by the OPAL detector in the years 1990 
through 1995 at e^e~ centre-of-mass energies close to the Z° resonance. Only data for 
which the tracking chambers and the electromagnetic calorimeter were fully operational 
are retained. The topology of e"'"e~— >■ t~^t~ events is characterised by a pair of back- 
to-back, narrow jets with low particle multiplicity. Jets are defined in this analysis by 
grouping tracks and electromagnetic clusters into cones with a half- angle of 35°, where 
each cone is assumed to contain the decay products of one of the tau leptons. Tau- 
pair candidates are selected by requiring two low-multiplicity jets (A^track < 6) with an 
average polar angle of |cos6'ave| < 0.68. Background events from other lepton pairs, 
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hadronic Z° — > qq decays, and four-ferinion events are reduced with cuts on the event 
topology and total visible energy. These selections produce a sample of 110326 tau-pair 
candidates with background f'^°'^~'^ = 1.56 ±0.10%, estimated by data control samples 
and Monte Carlo background samples. The tau-pair selection procedure is described 



in detail in reference [TO. 



One-prong tau decays are selected from this sample by choosing tau decay cones 
containing only one well-reconstructed charged track, consistent with coming from the 
origin. Tracks in this sample are required to have at least 40 jet chamber hits used 
in the calculation of the measured dE/dx, to have at least 3 hits in the z— chambers, 
and to have reconstructed momentum between 2 and 50 GeV/c. A sample of 143528 
candidates is selected by this procedure. 

The number of charged kaons in the sample is determined by maximising a likeli- 
hood function based on dE/dx. The r^^ Uj. K^ > Oh° branching ratio is then obtained 
by correcting the number of kaons for backgrounds and selection inefficiencies. As a 
cross-check of the general analysis procedure, the r~— ;> Uj. it~ > Oh'' branching ratio is 
also determined in an analogous fashion, and compared to the world average. 



5 Inclusive r ^ z/^ K > Oh^ Branching Ratio 



The tracks in the one-prong tau decay candidate sample consist of electrons, muons, 
pions, and kaons from both tau decays and non-tau sources. To determine the number 
of charged kaons in the sample, a maximum likelihood fit to the ionisation energy loss 
distribution of the tracks in the sample is performed. The likelihood function used in 
the fit is: 



C = exp 






j=l,N k=e,fj,,n,K 



(2) 



where W^^ is the dE/dx weight of charged particle j under particle hypothesis k, 

2" 



W''^ 



and where 



1 



27t:s{Pj)SresO-kj 



exp 



2 I s{Pj)SresO-kj 



(3) 



fk is the fraction of particle type k in the sample, where k is either pion, kaon, muon, 
or electron. The values of fk are constrained to be non-negative, and their sum 
is constrained to be 1. 



/' and A/' are the fraction, and uncertainty on the fraction, of muons in the sample, 
as estimated by Monte Carlo generated events. The estimated muon fraction is 
corrected using information from data and Monte Carlo muon control samples 
that are selected using information from the OPAL muon chambers. This con- 
straint is added to the fit to help distinguish between pions and muons in the 
sample, which have quite similar dE/dx. 

N is the total number of particles in the sample. 

-^meas i^ the mcasurcd dE/dx of the j^^ charged particle. 

Dpj^ed is the predicted dE/dx for charged particle j, calculated with the OPAL param- 
eterisation under particle hypothesis k, as obtained from the measured dE/dx of 
charged particles in hadronic Z° -^ qq events. 

(Tfcj is the dE/dx uncertainty, calculated using the OPAL parameterisation, as ob- 
tained from the measured dE/dx of hadronic Z° -^ qq events. 

Sres is the multiplicative correction to a^j, as determined from the one-prong control 
samples as described in Section 2.1. 

s(/5) is the (3 dependent multiplicative correction to -Dpred, as determined from the 
one-prong control samples. 

/(0) is the (p dependent correction to the measured dE/dx, as obtained from the 
one-prong control samples. 

Efforts are made to obtain a dE/dx parameterisation for the tau decay environment 
that is optimal for many particle species over a wide range of momenta. However, it is 
possible that the dE/dx corrections described in Section 2.1 may be somewhat more (or 
less) optimal for pions than they are for kaons in the momentum range of interest. Thus, 
to correct for any possible species-dependent quality differences in the parameterisation 
of dE/dx, an extra factor, Ck, is allowed to multiply the kaon predicted energy loss, 
-^^red) and is allowed to vary freely in the fit using the likelihood function found in 
Equation |. The value of Ck returned by the fit is Ck = 0.9943 ± 0.0009. 

Independent likelihood fits are performed in 13 momentum bins of variable size in 
the range of 2 to 50 GeV/c. A test of the fit with Monte Carlo generated samples 
verifies that the resulting estimates for the kaon fraction are efficient and unbiased. 

The number of charged kaons and pions found in the sample, summed over all 
momentum bins, is given in Table |I|, and the stretch dE/dx distribution of tracks in all 
momentum bins of the sample is shown in Figure p. Stretch energy loss. Si, is defined 
as: 

Si = (i^mcas " -^prod)/'^) 



where i is the particle hypothesis used to calculate the predicted dE/dx. The nor- 
mahsation of the predicted distributions of the kaons, pions, muons, and electrons in 
Figure ^is obtained from the results of the likelihood fit. The momentum distributions 
of the charged kaons and pions in the sample, as estimated by the likelihood fit, are 
shown in Figure |^. 



5.1 The Branching Ratio Calculation 

To determine the r~— ^i/^-K^ > Oh° branching ratio, the number of charged kaons 
found by the dE/dx likelihood fit in the one-prong candidate sample is corrected for 
backgrounds and selection inefficiencies: 

Bk-,„„.. = Br(. ^ ., K > Oh ) = ,,^_^„^„^^(i^^„„..-.) . (4) 

where 



A^^ is the number of tau decay candidates. There were 220652 tau decay candidates 
recorded in the barrel region of the OPAL detector between the years 1990 and 
1995. 

j-non-T jg ^j^g estimated background from non-r sources in the tau decay candidates 

(^ynon-r ^ 1.56 ±0.10%). 

A^TOTAL is the total number of kaons in the one-prong tau decay sample, as estimated 
by the likelihood fit. The number of kaons summed over all momentum bins in 
the sample is shown in Table |I|. 

-^bkgnd is *^^ number of background kaons in the one-prong tau decay sample, as 
estimated by Monte Carlo generated tau decay and hadronic Z° -^ qq samples. 
The number of background kaons summed over all momentum bins in the sample 
is shown in Table |I[ 

eK->oho is the efficiency for t'^Uj. K^ > Oh° decays in the tau-pair candidate sample 
to contribute kaons to the one-prong tau decay sample, and is shown in Table |I[ 
The efficiency has been corrected for biases introduced by the tau-pair selection 
procedure, and the uncertainty on the efficiency includes the systematic uncer- 
tainty arising from this correction. The bias corrections are all consistent with 1 
to within one standard deviation of their Monte Carlo statistical uncertainties. 

The result of the branching ratio calculation is: 

BK->ohO = 1.528 ±0.039 ±0.040%, 

9 



where the first uncertainty is statistical and the second is systematic. A summary of 
the systematic uncertainties, all of which will be described in subsequent sections, is 
shown in Table ^ 

As a cross-check of the general analysis procedure, we also determine the 
T'^UrTi^ > Oh° branching ratio. Thus, in a similar manner to the r~— ^z/^-K" > Oh'' 
branching ratio calculation, the total number of charged pions found in the one- 
prong tau decay sample by the likelihood fit is corrected for backgrounds and 
the T'^Urii^ > Oh° efficiency to yield B7,->ohO = 48.26 ± 0.25%, where the un- 
certainty is statistical only. The linear correlation coefficient between the OPAL 
T^^ Ur K^ > Oh° and r~^ Ur vr^ > Oh° branching ratios is approximately —20%. The 
OPAL t^^Ut-71^ > Oh'' branching ratio is consistent with the world average, derived 
from the branching ratios of the tau to the various exclusive final states that include a 
charged pion listed in reference Jl^, B^->ohO = 48.36 ± 0.23%, and thus we conclude 
that biases introduced by the analysis procedure are negligible. 



5.2 dE/dx Systematic Studies 

A significant source of systematic uncertainty in the measurement of the inclusive 
branching ratio is the uncertainty in the parameterisation of the predicted energy loss. 
To assess this systematic, we determine the sensitivity of the likelihood estimates of the 
number of kaons within the data sample to the uncertainties in the dE/dx correction 
factors obtained from the one-prong control samples. The likelihood function shown 
in Equation S is therefore modified such that: 



£' = exp (-i (n' - n)^V-^(n' - n)) C, 



where 



n' is a vector containing the central values of the dE/dx correction factors that describe 
the functions /(0), s(/5), and Sj-es, where the central values are as obtained from 
the fits to the one-prong control samples. 

n is a vector containing the dE/dx correction factors used in the likelihood fit. 

V is the covariance matrix for the dE/dx correction factors, as obtained from the fits 
to the one-prong control samples. 

A 'prior-belief method is used to assess the systematic uncertainty; in the first 
iteration of the fit, the correction factors are allowed to vary, and the returned values 
are found to be consistent with the values obtained from the control samples. In the 
second iteration, the likelihood fit is repeated, this time keeping the dE/dx correction 
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factors fixed to the values from the first iteration. The systematic uncertainty in 
A^TOTAL (and A^total) due to the parameterisation of dE /dx is then obtained from 
the square root of the quadrature difference of the fit uncertainties in A^total (^total) 
from the two iterations, and is shown in Table 0. 



5.3 Efficiency Estimation 

The efficiencies for kaons and pions from signal events in the tau-pair candidate 
sample to enter the data sample are estimated using signal events generated with the 
KORALZ 4.0 Monte Carlo generator and the TAUOLA 2.4 decay package, as described 
in Section 3. The efficiency for r"^ Vj. K~ > Oh'' decays in the tau-pair candidate 
sample to contribute to the one-prong sample, eK->oho, is calculated via: 

eK->oho = /^K- + (1 — /)eK->iho, 

where / is the ratio of the exclusive r^^ Vr K^ branching ratio to the inclusive 
r^^z/^K^ > Oh° branching ratio, as obtained from OPAL data (see Section 6.2). The 
eK">iho efficiency is calculated from the average of the relevant efficiencies appearing in 
Table ||, weighted by the associated world average branching ratios. The inclusive and 
exclusive branching ratios determined by this analysis are all less than 15% correlated 
with the world average branching ratios used in this efficiency calculation. All effi- 
ciencies have been corrected for biases introduced by the tau-pair selection procedure, 
and, in addition, the Monte Carlo efficiency estimates for the signal events to pass the 
jet chamber and z-chamber hit requirements are corrected using a control sample of 
data one-prong tau decays. The uncertainty on each efficiency includes the systematic 
uncertainties arising from these corrections. 

The efficiency for r^^ Vr tt^ > Oh° decays in the tau-pair candidate sample to con- 
tribute to the one-prong sample, e^->ohO, is calculated using the same procedure, ex- 
cept that the ratio of the r~-H>z/^7r~ branching ratio relative to the r~— s-i/t-tt" > Oh'' 
branching ratio is taken from world averages. 

The Monte Carlo modelling of the probability for photons to convert to an e"'"e~ 
pair in the OPAL detector material has been found by previous studies of data control 



samples to be correct to within an absolute uncertainty of about ±0.5% |18]. To assess 
the systematic uncertainty in the branching fraction associated with this modelling, we 
determine the change in the branching fraction when the efficiency for each exclusive 
final state is changed by +0.005 for each associated photon in the final state; thus the 
efficiency for r~— >i/,- K~7r° decays is changed by 2 x 0.005 = 0.01, and the efficiency for 
r~— s-i/^K^tt^tt" decays is changed by 4 x 0.005 = 0.02. We also determine the change 
in the branching fraction when the efficiency for each exclusive final state is changed 
by —0.005 for each associated photon in the final state. The resulting systematic 
uncertainties in the branching fractions are shown in Table 0. 
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The variation of the branching ratio due to alternative intermediate resonant 
structure scenarios for the t^^Ut-K^tt^tc^ final state is assessed using efficiency 
estimates from signal events generated by the modified version of TAUOLA 2.4 
through non-resonant production only. The efficiency calculated from this sample for 
r^^ Ur K^7r°7r° to contribute to the data candidates is within OAa of the central value, 
where a represents the combined Monte Carlo statistical uncertainties. Since the two 
different scenarios produce consistent results for the efficiency, and since there is no 
a priori reason to expect a strong efficiency dependence on the intermediate resonant 
structure, no systematic uncertainty for this effect is assigned. 

In a similar fashion, the variation of the branching ratio due to alternative inter- 
mediate resonant structure scenarios for the r~— s^z/,- K~K°7r° final state is explored 
using events generated through non-resonant production only. The r~— >i/,-K~K''7r° 
efficiency calculated from these events is within 0.6a of the central value, thus no 
additional systematic uncertainty is assigned. 

Since the observed momentum distributions for both the pions and the kaons in 
the sample are well described by the Monte Carlo generated events, no systematic 
uncertainty is assigned for possible error in either the estimate of the total efficiency, 
or the momentum dependence of the efficiency estimates, since it is apparent that such 
errors are small relative to the other systematic uncertainties, as listed in Table 0. 



5.4 Background Estimation 

The number of background kaons in A^total is estimated by applying the same 
selection criteria to Monte Carlo samples of hadronic Z^ -^ qq decays and three-prong 
tau decays. From the number of selected events that contain kaons, the estimated 
number of background kaons is derived, as summarised in Table |l|. There are very 
few background kaons in A^total^ most of which result from the small amount of 
low-multiplicity hadronic Z^ -^ qq events contaminating the tau-pair sample. Due to 
this small contamination, the systematic uncertainty associated with the background 
correction is neglible. 

The number of background pions in A^total is estimated in a comparable manner. 



6 Exclusive Branching Ratios 

To distinguish between one-prong tau final states with and without accompanying neu- 
tral particles, information from the electromagnetic calorimeter is used. Tau decays 
that include neutral particles tend to have a larger relative electromagnetic energy 
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deposition, E/p, associated with the tau decay cone than tau decays without accom- 
panying neutrals. Such decays also tend to have a larger number of associated clusters, 
iVcius, in the electromagnetic calorimeter than tau decays without neutrals. To en- 
sure accurate modelling of the E /p and A^cius distributions of one-prong tau decays 
to charged kaons, corrections to the Monte Carlo simulation are obtained from a data 
T^^v^ 71^ > Oh° control sample. The following sections describe these corrections, and 
how they are applied in the determination of the r~— i^z/^K~ branching ratio. 



6.1 Corrections to the Modelling of E/p and N^us 

The two-dimensional (p, E/p) and (p, A^cius) distributions of data r^^ z/^ vr^ > Oh° de- 
cays are used to determine corrections to the Monte Carlo modelling of the E/p and 
A^cius distributions of t~^ z/^ K^ > Oh° decays. To this end, the data one-prong tau 
decay sample is divided into 13 bins of momentum between 2 to 50 GeV/c and 6 bins 
of E/p, and the number of charged pions and kaons within each bin is determined 
by maximising the dE/dx likelihood function shown in Equation ^. After background 
correction, this procedure yields the (p, E/p) distributions of r^^ Ur vr^ > Oh° and 
t^^Ut-K^ > Oh° decays. The data one-prong sample is also divided into 13 bins of 
momentum and 4 bins of Ndns, smd the [p, Ndns) distributions of r~^ z/^ vr" > Oh'^ and 
r~— s-z/^K" > Oh° decays are determined in a comparable manner. 

The momentum distribution of the charged pions in the one-prong tau decay sam- 
ple, as shown in Figure |^, appears to be adequately modelled by Monte Carlo generated 
r~^z/^7r~ > Oh° events. However, the E/p and A'^dus distributions, obtained by sum- 
ming the {p, E/p) and {p, Ndns) distributions over momentum bins, are not as well 
modelled by the Monte Carlo, as is shown in Figure ^. 



6.1.1 Corrections to E/p 



Separate corrections are derived for the Monte Carlo E/p modelling of r^— i>z/^7r^, 
r~— >■ Ur 7r~7r°, and t~^ z/^ 7r~7r°7r° decays. These three decay modes produce 95% of the 
charged pions found in the one-prong tau decay sample. The E/p corrections obtained 
for t^^Ut-71^ are applied to the Monte Carlo modelling of the E/p of r^^z/^-K", 
and, in a similar fashion, the corrections to the E/p distribution of r~^z/^7r^7r° 
(r~— >z/T-7r~7r°7r°) decays are applied to that of r~-^z/T-K~h° {T~—>-UrK~ > 2h'^) de- 
cays. 

The corrections to the r~— ^z/t-tt", r~-^z/T-7r~7r°, and r~^z/T-7r~7r°7r'' Monte Carlo 
modelling of E/p are obtained by first correcting the number of charged pions within 
each {p, E/p) bin of the one-prong tau decay sample for background pions from 
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hadronic Z° -^ qq events. The number of pions in each bin is also corrected for the num- 
ber of pions from tau sources other than r^— i> Uj. vr", r^— i> z/^- tt^tt^, and r^^ z/^ tt^ti^ti^ 
decays using world average branching ratios and information from Monte Carlo gener- 
ated events. The following x^ function is then minimised: 



x'=j:j: 






(5) 



where the index i runs over the momentum bins, the index j runs over the E/p bins, 
and the index / runs over t^^Ut-tt^ , r^^z/,- 7r~7r°, and r^— ^z/^vr^vr'^vr^. Also: 

N[j is the number of charged pions, after corrections for non-r^^z/^ tt^ < 27r*^ back- 
ground, found by the dE/dx likelihood fit in the {i,jy^ bin of the {p,E/p) 
distribution of one-prong tau decays. 

AN^j is the statistical uncertainty on N^j, as obtained from the dE/dx hkelihood fit. 

Bi is the world average branching ratio of tau decay /, taken from reference |jl^ . 



ei is the efficiency for events from tau decay / in the tau-pair candidate sample to 
contribute to the one-prong tau decay sample. 

Slj is the probability, calculated using Monte Carlo generated events, for tau decays 
of type / in the one-prong tau decay sample to contribute to the (i, j)*^ bin of 
the (p, E/p) distribution. 

Sj is the correction to the j^^ bin of the Monte Carlo E/p distribution of tau decays 
of type /. The iSj are the only parameters allowed to float freely in the flt. 

s\ is a normalisation factor used to ensure that the correction to E/p does not affect 
the shape of the momentum distribution of decays of type /, and is determined 
from the condition J2j s[SljSj = P/, where P- is the probability, as calculated 
from Monte Carlo generated events, for tau decays of type / to contribute to the 
i^^ bin of the momentum distribution. 



6.1.2 Corrections to A^cius 

Corrections, T, to the Monte Carlo modelling of A'^^us of r~— > z/^ vr", r~— > z/^ 7r~7r'', and 
r~— i> z/t- 7r~7r°7r° decays are obtained in a comparable manner to the method used to ob- 
tain the E/p corrections. Again, the r^^ z/^ vr^, r^— > i/^ 7r~7r°, and r~— > z/^ tt'tt'^tt^ N^us 
distributions are similar to those of T^^UrK^, r^^z^T-K~h°, and r^—'-UrK^ > 2h° 
decays, respectively. 
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The absolute linear correlation coefficients between A'^dus and p, and A^cius and E/p 
are less than 20% in all classes of signal. Thus, the A^'^us and the (p, E/p) distributions 
are used as discriminators in this analysis, and A'cius is treated as a variable uncorrelated 
to both p and E/p. 



6.2 The r ^z/^K Branching Ratio 

The r~— >z/t-K~ branching ratio, Bk-, is determined by a x^ fit to the {p,E/P) and 
A'^cius distributions of one-prong tau decays with a charged kaon, where the modelling 
of the distributions is corrected using the results of the previous section. The following 
X^ is minimised: 



X 



+ 



EE 



m-N^^'^-EieisiSl^G^R 



tj J 






E 

k 



Nf 



Nt^'-'-TieiTiHiRf 



AiVf 



(6) 



where / is an index referring to either r~— >z/t-K^, r^— >z/T-K^h°, or r^^z/^-K^ > 2h'', 
and the indices i, j, and k run over the p, E/p, and A'^dus bins, respectively. The 
number of background kaons within each bin is A^^'^snd^ while AA^^ is the statistical 
uncertainty on the number of kaons within each bin, A^^. Also: 
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/) 



Bx-ho + Bk- 
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and 



/ is the ratio of the OPAL r ^ z/^ K branching ratio to the OPAL r 
inclusive branching ratio, and is allowed to float in the fit. 



•V. 



V K- > Oh° 
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-^TOTAL is the total number of charged kaons in the one-prong tau decay sample, as 
determined in Section 5, and shown in Table 0. 

e; is the efficiency for events from tau decay / in the tau-pair candidate sample to 
contribute to the one-prong tau decay sample. The ex-ho, eK->2hO; and eK->ihO 
efficiencies are calculated from the average of the relevant efficiencies appear- 
ing in Table ^, weighted by the associated world average branching ratios. All 
efficiencies have been corrected for biases introduced by the tau-pair selection 
procedure. 

Bx-ho and BK->2ho are the world average branching ratios for r~— >z/,-K~h° and 
r~— »• Ur K~ > 2h*^ decays, respectively, derived from results presented in refer- 
ence Wi 



S\j is the probability, calculated using Monte Carlo generated events, for tau decays of 
type / in the one-prong tau sample to contribute to the {i^jY^ bin of the (p, E /p) 
distribution. 

S'j , Sj '^ , and SJ '^ '" are the corrections to the j^^ bin of the Monte Carlo E /p 
distribution of the various one-prong tau decays to pions, as determined from the 
X^ fit using Equation ^ 

s\ is a normalisation factor used to ensure that the correction to E /p does not affect 
the shape of the momentum distribution of decays of type /. 



T^ is the probability, calculated using Monte Carlo generated events, for tau decays 
of type / in the one-prong tau sample to contribute to the /c**^ bin of the A^cius 
distribution. 



q-TT ^ j-TT -K ^ g^j^j -j-TT TT TT ^^^ ^]-^g corrcctlous to thc fc*^ bin of the Monte Carlo A^cius 
distribution of the various one-prong tau decays to pions, as determined in the 
previous section. 

The value of / derived from the result of the fit using Equation ^ is 
/ = 0.436 ± 0.013, where the uncertainty is statistical only. The t~—*Vj. K~ branching 
ratio derived from this result is: 

Bk- =Br(r-^z/^K-) = 0.658 ± 0.024 ± 0.029%, 

where the first uncertainty is statistical and the second is systematic, and the exclusive 
branching ratio is 60% correlated with the OPAL inclusive r~— s-i^t- K~ > Oh*^ branching 
ratio. A summary of the systematic uncertainties on Bk- is shown in Table ^. 

The statistical uncertainty in Bk- due to the statistical uncertainty in the number 
of kaons within each bin of the (p, E/p) and Nc\ns distributions is assessed by repeating 
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the fit many times, each time varying the central value of the number of kaons within a 
particular bin, a bin at a time, by plus, then minus, one statistical standard deviation. 
The square root of half of the quadrature sum of the net variations produced in the 
branching ratio is taken as the statistical uncertainty. 

The total x^ per degree of freedom of the fit using Equation ^ is 63.0/50, where 
56.6/46 results from the x^ per degree of freedom between the data and predicted 
(p, E /p) distributions. Figure |^ shows the E /p and A^cius distributions of one-prong 
tau decays with a charged kaon, where the normalisation of the contributions to the 
predicted distribution comes from the OPAL r^^ v^ K~ and r^^ Vj. K^ > Oh*' branch- 
ing ratios. If the E /p and A^cius corrections are not applied, the total x^ per degree of 
freedom from the fit is 79.9/50. 

To check for undue variation that may be produced by the choice of binning used 
in the fit, the above procedure is repeated using various binning schemes. The RMS 
variation of the results returned by the different fits is less than the statistical uncer- 
tainty from the original fit. In addition, a test of the entire fitting procedure (including 
the calorimeter variable correction procedure described in Sections 6.1 and 6.2) is per- 
formed with Monte Carlo generated samples, and verifies that the resulting estimates 
of / are efficient and unbiased. 



6.3 dE/dx Systematic Error 

The uncertainty in the N^ within each A^cius or (p? E/p) bin due to uncertainty in 
the dE/dx correction factors is determined as in Section 5.2. The dE/dx uncertainties 
are highly correlated between the bins. Thus the uncertainty in Bk- due to the dE/dx 
systematic uncertainty in the number of kaons within each bin of the (p, E/p) and 
A^cius distributions is assessed by repeating the fit many times, each time varying the 
the central value of the number of kaons within a particular bin, a bin at a time, by 
plus, then minus, one standard deviation of the dE /dx systematic uncertainty. Half of 
the linear sum of the absolute variations produced in the branching ratio is taken as 
the dE/dx systematic uncertainty in Bk- , and is shown in Table 13. 



6.4 Systematic Uncertainty Associated with the Efficiency Es- 
timation 

The uncertainty in the r~— ;> v^ K~ branching ratio due to the limited Monte Carlo 
statistics used to determine the eK-, CK-ho, and eK->2ho efficiencies is assessed via the 
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prior-belief method by modifying Equation ^ such that: 

where 

e'f. and Ae'^ are the efficiency and associated uncertainty, due to hmited Monte Carlo 
statistics, for kaons from signal channel k in the tau-pair candidate sample to 
enter the data sample. 

ejfc is the efficiency used in the fit for kaons from signal channel k to contribute to the 
data sample. 

In the first iteration of the fit the ek are allowed to fioat. In the second, they are 
fixed to the values from the first iteration. The square root of the quadrature difference 
in the errors in the branching ratio is taken as the systematic uncertainty associated 
with the limited Monte Carlo statistics used in the assessment of the efficiencies, and 
is shown in Table §. 

Various world average exclusive branching ratios are used in the calculation of ex-ho 
and eK->2hO (for instance, ex-ho is the average of ex-TrO and ex-x"; weighted by the 
world average branching ratios Bx-tt" and Bx-xo)- In addition, these same branching 
ratios are used to assess the inclusive r~— ^^'rK~h° and r^^z/^K~ > 2h° branching 
ratios appearing in Equation ^. The uncertainty in the r^^ v^ K~ branching ratio 
due to the uncertainties in these world average branching ratios is assessed by varying 
the branching ratios, one at a time, by plus, then minus, one standard deviation. 
The square root of half of the quadrature sum of the net variations produced in the 
branching ratio is taken as the systematic uncertainty, and is shown in Table ^ 

The systematic uncertainty in the efficiency and branching fraction due to the 
uncertainty associated with the Monte Carlo modelling of the probability for photons 
to convert to an e"^e~ pair in the OPAL detector material is assessed as described in 
Section 5.3. 



6.5 Systematic Uncertainty Associated with the ModeUing of 

E/p and iVcius 

The uncertainty in the t~^>-Vj. K~ branching ratio due to the statistical uncertain- 
ties in the E/p correction factors is assessed using the prior-belief method by modifying 
Equation ^ such that: 

x'' = x' + (n'-nfV-i(n'-n), 



where 



n' is a vector containing the central values of the E/p corrections, as obtained from 
the fit to the r~— s^z/t-tt^ < 27r° sample. 

n is a vector containing the values of the E/p corrections used in the x^ fit. 

V is the covariance matrix for the S correction factors, as obtained from the fit to the 
r~— >i/^7r^ < 27r° {p,E/p) and A^cius distributions. 

In the first iteration, all bins of the three sets of S corrections are allowed to vary 
in the fit, and the returned values are found to be consistent with the input values. 
In the second iteration, the fit is repeated, this time keeping the E/p corrections fixed 
to the values from the first iteration. The systematic uncertainty in Bk- due to the 
modelling of E/p is then obtained from the square root of the quadrature difference of 
the fit uncertainties from the two iterations. 

The uncertainty in the r~— i>z/^K~ branching ratio due to the statistical uncertain- 
ties in the A'^^us correction factors is assessed in an analogous fashion. The combined 
systematic uncertainty in Bk- due to the corrections to the modelling of E/p and A^cius 
is shown in Table 0. 



7 Summary and Discussion 

From an analysis of the ionisation energy loss of charged particles selected from a 
sample of 220652 tau decay candidates recorded in the barrel region of the OPAL 
detector at e"'"e~ centre-of-mass energies close to the Z° resonance, we determine the 
branching ratios: 

Br(r-->z/^K- > OhO) = 1.528 ± 0.039 ± 0.040% 

Br(r"^z/^K-) = 0.658 ± 0.024 ± 0.029% 

Br(r-^ z/^ K- > lh°) = 0.869 ± 0.031 ± 0.034% 

where the first uncertainties are statistical and the second are systematic. The 
t'^Ut-K'^ and r^^z/^K~ > Ih*^ branching ratios are 60% and 80% correlated to the 
r^— i> Ur K^ > Oh*^ branching ratio, respectively. 

Both the inclusive and exclusive branching ratios are in agreement with previous 
measurements, shown in Table |^, and both contribute significantly to the reduction of 
the uncertainty on the world averages of the respective quantities. 
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The theoretical prediction for the r ^ z/,- K branching ratio is calculated via [^ : 

Br(r-^..K-) = ^^/||V„.p (l - m|/m^)^ 

= 0.705 ±0.008%, (7) 

where the values of /k = 0.1598 ± 0.0015 GeV and |V„,| = 0.2196 ± 0.0023 are both 
taken from reference [0. The OPAL measurement is in agreement with this prediction 
to within 1.3 standard deviations of the combined uncertainties. 

The OPAL measurement of the t^—^Ut-K^ branching fraction may be used to de- 
termine either /k or \Vus\ through Equation |^. Assuming /k = 0.1598 ± 0.0015 GeV, 
a value of {Vusl = 0.2121 ± 0.0063 is obtained. Assuming a value of {Vusl = 
0.2196 ± 0.0023, a value of /k = 0.1544 ± 0.0046 GeV is obtained. In both cases, 
the uncertainty on the determined quantity is dominated by the uncertainty on the 
r~— >i/^K~ branching fraction. 

The theoretical prediction for the r~^z/T-K~ > Oh° branching ratio is obtained by 
summing the theoretical predictions for the five states appearing in Table ^: 

Br(r-^ Ur K- > Oh°) = 1.571 ± 0.043%. 

The OPAL measurement is in agreement with this prediction to within 0.7 standard 
deviation of the combined uncertainties. 

The OPAL collaboration has previously presented results for the tau branching 
fractions to the i^T-K~7r~7r+ > Ovr*^, Ut-tt'K^, and i/^7r~K° > Ivr^ final states [26, 27]. 



These branching ratios can be combined with the OPAL measurement of the r — > z/^ K" 



branching ratio, using isospin relations appearing in |]28[, to yield an estimate of the 



branching fraction of the r to final states with strangeness —1: 

Bt{t- -^ [strangeness=-l]) = 2.81 ± 0.16 ± 0.10%, 

where the calculation of the uncertainties takes into account all correlations between the 
various OPAL results. The calculation neglects contributions from r~-^ u^ K~ri to the 
final state, and assumes that the the r~— >z/T-K^7r^7r+ > Ovr^ final state is dominated 
by t^—i-i't-K'tt^tt^, and that the t" —>■ z/^tt^K*^ > Ivr*^ final state is dominated by 
T^ — i> UrTc^K^Ti^. The OPAL result is in agreement with the previously published 



result p9[, Br(r — > [strangeness=-l]) = 2.87 ±0.12%, to within 0.3o" of the combined 



uncertainties. 
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Nr 
xnon— T 

# in one-prong sample 


220652 

1.56 ±0.10% 

143528 


^^ TOTAL 
^^ bkgnd 
^■rr->Oh° 


75566 ± 389 ± 246 

89 ±16 

0.720 ±0.001 ±0.001 


^^ TOTAL 

^bkgnd 

eK->OhO 


2526 ± 64 ± 52 

9±4 

0.759 ±0.006 ±0.007 



Table 1: Pion and kaon composition of the one-prong tau decay candidate sample, as 
estimated by the likelihood fit to the measured dE/dx of tracks in the sample. The 
first uncertainty on iV^oTAL ^^^ ^total i^ the statistical uncertainty from the fit, 
and the second is due to the systematic uncertainty in the dE/dx correction factors. 
Also shown are the estimated backgrounds in A^tqtal ^^^ ^total; along with the 
efficiencies for signal events to contribute to the one-prong tau decay sample. The first 
uncertainty on each efficiency estimate is due to the limited statistics of the Monte 
Carlo generated samples, while the second is due to the uncertainty in the various 
branching ratios used to calculate the efficiency, as listed in Table 0. 



T 


Decay 


BR (%) 


Efficiency 


Mode 






T~ 


-^I/^TT 


11.08 ±0.13 


0.810 ±0.003 


T^ 


-^z/^7r"7r° 


25.32 ±0.15 


0.725 ±0.002 


T^ 


-^z/,7r-K0 


0.83 ±0.08 


0.604 ±0.007 


T^ 


-^z/^7r"7r°7r° 


9.15 ±0.15 


0.642 ± 0.003 


T^ 


-^Urii-K^K^ 


0.121 ±0.021 


0.454 ±0.016 


T^ 


-^z/.tt-KOttO 


0.39 ±0.05 


0.550 ±0.015 


T^ 


-^z/.tt-ttOttV 


1.11 ±0.14 


0.569 ±0.005 


T^ 


-^z/^7r"r7± > 07r° 


0.188 ±0.025 


0.655 ±0.024 


T^ 


-^z/^K- 


- 


0.843 ±0.010 


T^ 


-^z/^K^ttO 


0.52 ±0.05 


0.757 ±0.011 


T^ 


-^ Ur K-K° 


0.159 ±0.024 


0.661 ±0.013 


T^ 


-^ Ur K~r] 


0.027 ±0.006 


0.757 ±0.011 


T^ 


-^z/,K-7rV 


0.080 ±0.027 


0.664 ±0.025 


T^ 


-^z/,K-KV 


0.151 ±0.029 


0.544 ±0.012 


T^ 


-^UrK-n'n'n^ 


0.035 ±0.024 


0.643 ± 0.025 


r" 


-^z/.K-KVttO 


0.000 ±0.018 


0.560 ±0.012 



Table 2: Branching ratios and efficiencies used to estimate the efficiencies for the 
T^^UrK^ > Oh'^ and r^—^UrTi^ > Oh'' decays in the tau-pair sample to contribute to 
the one-prong sample. The r^—^Ur K'K^tt^tc^ branching fraction is estimated from the 
result appearing in |ll9| . 
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Branching Ratios (%) 


^7T->0h0 


BK->OhO 


Bk- 


Central value 


48.17 


1.528 


0.658 


a (stat) 


±0.25 


±0.039 


±0.024 


a {dE/dx sys) 


±0.16 


±0.032 


±0.017 


a {E/p and A^cius modelling) 


- 


- 


±0.014 


a (photon conversion modelling) 


±0.68 


±0.015 


±0.003 


a (MC Stat) 


±0.09 


±0.012 


±0.006 


a (efficiency sys) 


±0.05 


±0.014 


±0.017 



Table 3: Summary of the branching ratio central values and sources of uncertainties. 
The second-to-last uncertainty is the systematic uncertainty arising from the limited 
statistics of the Monte Carlo generated samples used to estimate the selection effi- 
ciencies. The last uncertainty is due to the systematic uncertainty in the efficiency 
correction arising from the uncertainties in the various world average branching ratios 



used in the weighted average calculation of the r 
and t^^Ut-K^ > lh° efficiencies. 



■u. 



K~ > Oh°, r 



■U^TT 



>0h" 



r 


- DECAY 
MODE 


THEORY 
BR (%) 


r~- 


-^z/^K^ 


0.705 ±0.008 13 


r~- 


-^z/.K-vrO 


0.48 ± 0.02 [0 


r~- 


^l^rK-K^ 


0.111 ±0.031 [0] 


r~- 


^iyrK-7T%^ 


0.12 ±0.02 [g 


T'- 


^z/,K-KV 


0.155 ±0.005 13 



Table 4: Theoretical predictions for tau branching ratios to the dominant one-prong 
final states that include a charged kaon. 
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T- DECAY 
MODE 



r 



■v^K- 



UrK- > Oh" 



EXPERIMENT 
BR (%) 



0.696 ±0.025 ±0.014 

0.85 ±0.18 

0.66 ±0.07 ±0.09 

0.59 ±0.18 

0.658 ±0.024 ±0.029 



ALEPH99I19 



DELPHI94|2| 



CLE094|4 



DELC084| 

(this analysis) 



1.52 ±0.04 ±0.04 

1.54 ±0.24 

1.60 ±0.12 ±0.19 

1.71 ±0.29 

1.528 ±0.039 ±0.040 



ALEPH99|T9[ 
DELPHI94|2| 



CLE094[|24| 
DELC084II 
(this analysis) 



Table 5: Summary of the r — >i/^K and r ^i/^K > Oh° branching ratios. 
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Figure 1: (a) shows the ionisation energy loss -Dprcd as a function of the momentum 
for various particles in the OPAL jet chamber, (b) shows the resolution power 7?.jj 
expressed in terms of the d-E/dx resolution o", for various pairs of particle species i and 

i- 
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Figure 2: (a) is the stretch dE/dx distribution under a pion hypothesis for tracks in the 
data T^—i-h'T- h^ > Oh° candidate sample (points). The overlaid curves are the predicted 
distributions for kaons, pions, muons, and electrons in the sample. The normalisation 
of the curves is obtained from the results of the likelihood fit described in Section 5. (b) 
is the distribution of the data points in plot (a) divided by the predicted distribution. 
The shaded area represents the approximate one sigma dE/dx systematic uncertainty 
envelope on the predicted distribution. The x^ per degree of freedom between the data 
and a line centred at y = 1, taking into account both statistical and dE/dx systematic 
uncertainties, is 36.5/47. 
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Figure 3: (a) and (b) are the momentum distributions of charged kaons and pions, 
respectively, in the data one-prong tau decay sample (points), as obtained from the 
results of the d£^/da; likelihood fit described in Section 5. The histograms are the 
momentum distributions predicted by Monte Carlo generated t^—>- Vj. K^ > Oh*^ and 
t~^Vt- Ti~ > Oh" events. The overall normalisation of the predicted distributions comes 
from the results of the dE/dx likelihood fit, while the relative normalisation of the 
exclusive decay modes contributing to the sample is taken from the relevant world 
average branching ratios appearing in reference [0. 
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Figure 4: (a) and (b) are the uncorrected E/p and Ndns distributions, respectively, 
of Monte Carlo generated t'^Ut-tt" > Oh" decays (histograms), along with the E/p 
and A^cius distributions of data one-prong tau decays with a charged pion in the final 
state (points). The relative normalisations of the various exclusive tau decays that 
contribute to the sample are taken from the relevant world average branching ratios 
listed in reference [^]. Backgrounds are on the order of 0.2% or less and are neglected 
in the plots. The last bin in each of the distributions contains overflow events. 
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Figure 5: (a) and (b) are the corrected E/p and A'cius distributions, respectively, of 
Monte Carlo generated r"^z/T-K~ > Oh° decays (histogram), along with the E/p and 
iVcius distributions of data one-prong tau decays with a charged kaon in the final state 
(points). The overall normalisation of Monte Carlo distributions is obtained from the 
results of the dE/dx likelihood fit using Equation §, while the normalisation of the 
T^^UrK^ component is obtained from the results of the x^ fit using Equation ^ 
Backgrounds are on the order of 0.4% or less and are neglected in this plot. The last 
bin in each of the distributions contains overflow events. 
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